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ABSTRACT 



The objective of this thesis was to simulate steady, uniform, 
laminar flow over a plane surface with zero pressure gradient. De- 
sign parameters and requirements for the flow field and construction 
of the free surface water channel are discussed. Data acquisition was 
accomplished by means of the Hydrogen Bubble Technique; and, once 
acquired was used to obtain a normalized velocity profile. An un- 
certainty analysis on the acquisition method and a comparison of the 
experimental results with the solution obtained by Blasius are 
presented. 



L 1 • JL 

NA/AL POSTGRADUATE SGJ!OOL 
MONTEREY, CALIF. 93940 





TABLE OF CONTENTS 




Section 




Page 


1. 


Introduction 


9 


2. 


Equipment 


13 


3. 


Procedure and Presentation of Data 


16 


4. 


Remarks and Recommendations 


20 




Bibliography 


22 




Figures 


24 



3 



LIST OF FIGURES 

Figure p age 

1. Similar Velocity Profiles in the Boundary Layer 24 

2. Blasius 1 Solution for the Nondimensionalized 

Velocity Profile in the Boundary Layer 25 

3. Photographs of Uniform Stream Profiles 26 

4. Flow Production Section 27 

5. Flow Production Section 28 

6. Test Section with Traversing Mechanism 29 

7. Test Section with HBT Apparatus 30 

8. Block Diagram of Hydrogen Bubble Generating 

Equipment 31 

9. Photographs of Generated Time Lines in the 

Boundary Layer 32 

10. Instantaneous Position of Bubble Lines in 

the Bubble Plane 33 

11. Comparison of Experimental Velocity Profile 

with the Blasius Solution 34 



5 



NOMENCLATURE 
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boundary layer thickness 

distance between time lines measured downstream 

functional variable for the nondimensionalized 
coordinate system 

Froude number 

acceleration due to gravity 
depth of water in the channel 
pulse wire pulse repetition rate 
pulse wire pulse duration 
nondimensionalized vertical distance 
Reynold's number along a plane surface 
Reynold’s number based on plate length 
scale factor 

average instantaneous downchannel velocity 

uniform stream velocity 

vertical component of velocity 

uncertainty in ZS. x 

uncertainty in m 

uncertainty in S 

uncertainty in u 

direction taken downchannel 

direction taken vertical to the flow and 
perpendicular to x 
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1 . 



Introduction. 



The boundary layer concept was introduced by Prandtl [1] in 1904. 
Four years later Blasius f 2 ] , using Prandtl*s boundary layer equations, 
obtained a solution for steady, uniform, laminar flow, with constant 
fluid properties and zero pressure gradient, over a plane surface. As 
cited in Schlichting [3], experimental verification of Blasius* solu- 
tion followed in 1924 by Burgers T4] and van der Hegge Zijnen [5], in 
1930 by M. Hansen [6], and most carefully by Nikuradse f 7 ] , in 1942. 

A brief review of the applicable boundary layer equations character- 
izing the classic Blasius problem is presented. The use of the assump- 
tions: 

(1) incompressibility, 

(2) constant fluid properties, 

(3) steady two-dimensional laminar flow, 

(4) zero pressure gradient, 

(5) and the usual boundary layer assumptions, 
results in reducing the full Navier-Stokes equations to the more 
tractable form (Schlichting [3]), 



u & + 



V 41) 



X 





O 



Employing the similarity hypothesis advanced by Blasius f2] (velocity 
profiles in the boundary layer are similar, see Figure 1) one has 
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Where 



u = - ft<\) 

/% 

<\ - 1 - ^Rx 

6 x 

and 

Rx = Ujt 

O'* 

Using the stream function for uniform flow and the aforementioned 
definitions Blasius obtained the following normalized differential 
equation of the boundary layer: 

f(»0P'(fp + t = O 

and solved it for the boundary conditions: 

(1) f’(oo ) = 1, 

(2) f (0) = 0, 

(3) f’(0) = 0. 

The results of Blasius 1 solution are illustrated graphically in Fig- 
ure 2 • 

Over the succeeding years, subsequent investigative permutations 
of Blasius 1 original problem were made. One such permutation, under 
accelerated scrutiny, was the study of heat transfer through a fluid 
stream undergoing a liquid to solid phase transformation. The freezing 
of the fluid in the flow field over a plane surface was examined mathe- 
matically by Lapadula and Mueller [8] and Beaubouef and Chapman [9]. 
Pertinent experimental information, for the flat plate geometry, was 
provided by Savino and Siegel [10,11]. It was noted that Siegel and 
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Savino omitted the effect of the flow field in the analysis of their 
problem. 

In order to augment the limited supply of experimental informa- 
tion and to include the effect of the flow field in the analysis, a 
thesis was undertaken by B. F. Nolan, at the Naval Postgraduate School. 
Nolan [12] used a flat plate similar in geometry but larger in size 
than was employed by his predecessors. The flow field effect was in- 
cluded in his analysis by assuming Blasius flow over the cooling plate 
which allowed the boundary layer to remain unaffected by the solid 
phase geometry. The establishment of the Blasius flow field for Nolan’s 
thesis was undertaken as an independent project. 

It was apparent that the two projects were interrelated the first 
generating construction and flow requirements for the second. The 
cooling plate assembly, insulation and instrumentation required that 
the bottom of a free surface water channel be the test section for the 
second proposal. This requirement eliminated the difficulties which 
would otherwise be generated by a leading edge, instrumentation and 
support assembly of a flat plate model suspended in the flow. The ice 
growth study further required that the boundary layer should be suf- 
ficiently large or the Reynolds number be rather small. Quantitative 
data acquisition has been and continues to be an inherent difficulty 
in low velocity water flow. Impact tubes have been impractical and hot 
wire or hot film methods required a relatively pure fluid medium, which 
was not the case in this instance. Hama and Nutant [13], in a study 
concerning the transition phenomenon, required low speed flows. They 
discarded the dye-injection method in favor of the hydrogen bubble techni- 
que (hereafter referred to as HBT) as a means of gaining qualitative and 
quantitative information about the flow field. 
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There have been a number of papers in recent years concerning the 
investigation of the HBT for flow visualization and measurement. Some 
of the more noteworthy of those reports were those of Clutter, Smith and 
Brazier [14] in 1959, Schraub, Kline, Henry, Runstadler and Littell [15] 
in 1965 and Hansen [16] in 1968. A perusal of those papers revealed 
that the HBT when properly employed, was capable of excellent qualita- 
tive and reasonable quantitative information. The reliability of the 
HBT with respect to data acquisition was a function of many variables, 
the more important of which were: 

(1) steadiness of the flow field, 

(2) generating wire diameter, 

(3) electronics, 

(4) photographic tenchique. 

It was the objective of this thesis to design a system which would 
generate the desired flow conditions and to obtain a velocity profile 
using the HBT. The experimental velocity profile was then to be com- 
pared with the analytic solution obtained by Blasius. 
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2. Equipment 0 



The apparatus was comprised of a free surface water channel, control 
tank, floor tank, pump and associated piping. Driving force for the flow 
was provided by a Byron Jackson Company Pump Model No.08H-74. The pump 
took suction at a large sunken floor tank and discharged through a 6 inch 
diameter pipe, containing a standard ASME flange type orifice meter, and 
thence through a transition piece (from a circular section of 6 inch 
radius to a 1 square foot cross-section) into the channel. The channel 
was divided into three parts, the flow production section, the test 
section and a connecting channel. After flowing through the channel, the 
water entered a large control tank where the energy of the flow was dis- 
sipated. The drain from the control tank was set to maintain a con- 
stant water level in the channel for a given flowrate. 

The flow production section had several features all of which con- 
tributed to the purpose of the system, that of generating a uniform 
laminar flow. As reported by Lumley f 1 7 ] , honeycomb of large length to 
diameter ratio could be used to effectively reduce background turbulence. 
The large length to diameter ratio allowed the flow to become fully 
developed in the cells thus to some extent damping longitudinal flow 
oscillations as well as those in the transverse and vertical directions. 
Two regions of honeycomb (L/d ratio 24 and 13.3, respectively) were used 
in this project and proved effective in eliminating vertical and cross 
channel fluctuations. Two wire screens were also placed downstream from 
the honeycomb material. 

The screens, placed in order of increasing mesh size, in conjunction 
with the honeycomb helped to damp out longitudinal oscillations in the 
stream and produced a uniform velocity profile (see Figures 3, 4 and 5). 
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In addition, an adjustable gate was situated at the head of the unit and 
two fine mesh screens were placed downstream of the test section. The 
latter features, in addition to a damping board, were added to minimize 
surface waves. 

In a study conducted by Binnie, Davies and Orkney [18], the produc- 
tion of a uniform stream was discussed in some detail. It was found that 
if care was taken to eliminate incoming disturbances and a Froude number 
(F = U/ V gh) of 0.5 was not exceeded, the velocity profile was indeed 
uniform. For our purposes a maximum Froude number of 0.02 was contem- 
plated as a design parameter. This value was not exceeded during the 
course of the data runs. Also, the desired uniform flows ranged from 
0.5 to 2.0 inches/second implying Reynolds numbers (R = U1/*V“ ) ranging 
from 4165 to 16,660, well below the critical Reynolds number for a flat 
plate (500,000) . 

The test section of the channel was 18 inches long, 11.5 inches wide 
and 12 inches high. The test surface had a rectangular array of holes 
(0.025 inches in diameter) tapped through it. These holes covered the 
area where the cooling plate of the associated project would be located. 
The array served as a means of insuring a constant head pressure over 
the test surface and also provided access for the HBT generating wire. 

A versatile traversing mechanism contained the HBT pulse wire electrode 
and provided a means of accurately locating the pulse wire perpendicular 
to the test surface. The test section including the traversing mechan- 
sim and HBT apparatus are shown in Figures 6 and 7. 

The report by E. Hansen [15] was consulted frequently during the 
design, construction and operation of the HBT equipment. A block dia- 
gram illustrating equipment arrangement and corresponding manufacturers 
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are shown in Figure 8. The triggered pulse amplifier, a near duplicate 
of the one discussed by E. Hansen, was constructed to provide pulse 
power to the generating wire. Fast response coupled with a wide range 
of voltage and current produced a good bubble line which was of prime 
consideration in the measurement of the velocity profile. 

The photographic equipment used to record the instantaneous bubble 
configuration in the test section was a Polaroid-Land Camera Model HOB, 
including a Polaroid Close-up Lens Kit no. 550, with Polaroid 3000 speed 
black and white film. A scale factor (S) was obtained by comparing a sub- 
merged scale in the bubble plane with that same scale out of water. With 
the bottom, top and opposite side of the test section shielded from ex- 
ternal light, the bottom was provided with an unshielded strip through 
which a collimated light source illuminated the bubble plane. 

The stream flowed into the control tank through a portion of the 
connecting channel. Regulation of the control tank and the water level 
was done with a parabolic gate valve which provided a linear change in 
the flowrate with a change in gate position. 
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3. 



Procedure and Presentation of Data. 



Once a steady flowrate and water level were established in the 
channel, necessary steps were taken to acquire a velocity profile. As 
stated previously, illumination of the bubble plane was provided by a 
collimated light source shining up through the test surface. The next 
step was to generate timed marker lines into the stream. 

The pulse repetition rate (m) and pulse duration (m,) were adjusted 

d 

so as to give optimum line definition and spacing. Photographic require- 
ments (f/4.57, shutter speed of 1/125 sec., close up lenses 4*2 and 44, 
range at 3 feet) were found prior to beginning the run. The use of the 
Polaroid-Land Camera allowed quick qualitative information on the adequacy 
of the photographs. Typical photographs are shown in Figure 9. 

Longitudinal spacings (Ax) between marker lines (for vertical 
channel positions ranging from 0.0 to 1.2 inches) in 0.1 in. intervals 
were measured. The velocities were calculated from the formula, u = 

Ax*m*S, where, as stated previously, m is the pulse-wire pulse repti- 
tion rate and S is the scale factor. To calculate the normalized velo- 
city profile we require the quantities u/U and £ , where 

SL - AJL 

U AX 

and 

- ys 

The boundary layer thickness ( & ) has been defined (Schlichting 3) 

as y = £> for u/U = 0.99. This corresponds to = 5.0 or £ = 

y/5.0. Thus, for each y, in the range indicated, normalized velocity 
and vertical distances were obtained through the measurement of the 
variables graphically defined in Figure 10. 
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Figure 10 is presented as a diagramatic representation of the in- 
stantaneous position of the marker lines in the bubble plane and should 
not be misconstrued as depicting velocity profiles. The inflection point 
seen in the marker lines near the wire was not representative of the flow 
field; but, was a characteristic phenomena exhibited by bubbles moving 
in the boundary layer. The presence of the inflection point may be ex- 
plained as follows: the size of the bubbles generated by the pulse wire 

were in the order of the wire ® s diameter. Thus, the bubbles that col- 
lected within the wake of the cylindrically shaped wire fell in a re- 
gion of the flow field where there was considerable momentum deficiency 
due to the presence of the wire. Since there is a vertical component 
of velocity within the boundary layer (Schlichting [3]) the bubbles 
were convected upwards and more or less parallel to the wire. An addi- 
tional contribution to this convective motion came from the secondary 
fluid motion generated within the wake. It is therefore important to 
note that one should not take data very near to the generating wire 
and should also be extremely careful in this and other studies regarding 
the stability of flows as far as the inflection points are concerned. 

It was remarkable indeed that such a small wire could, together with the 
aid of the bubbles, gave the false appearance of the existence of an 
inflection point when, in reality, there was not one. 

In following the guidelines set down by Kline and McClintock [19], 
and Kline, Schraub, Henry, Runstadler and Littell [14], it became ap- 
parent that the significant uncertainty in u (W ) was the uncertainty in 
Ax (W ^ x ) o Since the bubbles were of sufficiently small size and 
photographs were taken over seventy diameters downstream of the pulse 
wire, bubble rise and velocity defect were minimally contributing 
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uncertainties and therefore omitted from the analysis. Great care was 
taken to produce a uniform laminar flow field and to place the generating 
wire vertically in the stream and perpendicular to the test surface. 

The purpose of this was to minimize the uncertainties induced by a 
transverse slope in the bubble plane so that such errors could also be 
omitted from the analysis. 

The velocity in the bubble plane was given by 

U = ax • m • s 

It followed then, that 




Where W and W were the uncertainties in the pulse repetition rate 
ms 

and scale factor, respectively. As stated by Shraub, Kline, et. al. [15] 

W and W may be neglected since they were small relative to W 
ms A ^ 

Therefore , 

Uu - 

Normalizing, one has 

uk - 

u 

A typical x for the runs made with this equipment was 0.30 + 0.033 
inches (20:1). Therefore, 

l/)u - IjOax- ~ o. 033 O.U 

O AX 0.3o 



-\ f ( • S - vJa*) — V>^ax 
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AX 



AX 
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or Che percent uncertainty in u was 11%. 

A comparison of the normalized velocity profile, reduced from the 
experimental data, with the Blasius solution shows a maximum deviation 
of 12%. Hence, it seemed reasonable to conclude that the accuracy of 
the apparatus at present was in accord with the calculated uncertainty. 

The normalized velocity profile, obtained experimentally, was 
shown in Figure 11, for three different flow rates. The experimental 
profile appeared to be in agreement with that theoretically obtained by 
Blasius indicating that the desired flow field had been established. 
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4. Remarks and Recommendations. 

A more accurate velocity profile could be obtained if the un- 
certainty in A x were reduced. A sizeable reduction in this uncertainty 
could be achieved if there were better bubble line definition. The most 
aggravating contribution to the lack of a better definition is the build 
up of impurities on the generating wire. This build-up generates stream 
ing sites on the wire and limits its useful age. As a result of the 
streaming phenomenon we get a thicker time line and consequently a 
larger uncertainty interval. 

The following recommendations are submitted in order to improve 
marker line quality and generally enhance the efficiency of the equip- 
ment: 

(1) Provide the triggered pulse amplifier with a polarity re- 
versing switch. A series of pulses of reversed polarity 
to the generating wire has a cleansing effect. 

(2) Provide a filtering system for the water in the floor tank. 

This would decrease the number of impurities and may even make 
practical the use of a hot film device as a second means of 
measuring velocity. 

(3) Change from Nichrome alloy wires to high tensile strength 
Platinum alloy wires as the Platinum wires have a longer 
useful lifetime and tend to resist the buildup of impurities. 

(4) Replace the present power amplifier and pulse generator ar- 
rangement with a Hewelett-Packard pulse generator Model No. 
214A, with 200 watt output. This pulse generator coupled with 
the triggered pulse amplifier should give a more versatile 
pulse characteristic. 
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(5) Construct an adjustable constant head tank ahead of the flow 
production section to eliminate any pump surging. 

(6) Provide a calibrated V-notch weir at the end of the connecting 
channel as a second means of measuring flowrate. This would 
also tend to minimize surface waves. 

(7) Try the perforated plate and rubberized hair suggested by 
Nutant [12], in the flow production section. 

(8) See if it is feasible to apply a correction factor for bubble 
motion in the boundary layer. The fact that the normalized 
velocity profile is in excellent agreement outside the bound- 
ary layer and falls off markedly on entry into the boundary 
layer deems a correction factor appropriate. 
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FIGURE 1. SIMILAR VELOCITY PROFILES IN THE BOUNDARY LAYER 
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U = 0.81 INCHES/SECOND , W - 1.85 GPS 



FIGURE 3 . PHOTOGRAPHS OF UNIFORM STREAM PROFILES 
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FIGURE $. FLOW PRODUCTION SECTION 



TRAVERSING MECHANISM 
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FIGURE 6. TEST SECTION WITH TRAVERSING MECHANISM 




PULSE GENERATOR 
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FIGURE 7. TEST SECTION WITH HBT APPARATUS 
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FIGURE 8. BLOCK DIAGRAM OF HYDROGEN BUBBLE GENERATING EQUIPMENT 












u = 1.05 inches/ second, 'm= 2,78 cps 
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FIGURE 9, PHOTOGRAPHS OF GENERATED TIME LINES 
IN THE BOUNDARY LAYER 
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FIGURE 10. INSTANTANEOUS POSITION OF BUBBLE LINES IN THE BUBBLE PLANE 
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